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Abstract 



I examine the determination of the QCD color factor ratio Ca/Cf from the scale 
evolution of particle multiplicity in e + e~ three jet events. I fit an analytic expres- 
sion for the multiplicity in three jet events to event samples generated with QCD 
multihadronic event generators. I demonstrate that a one parameter fit of Ca/Cf 
yields the expected result Ca/Cf=2.25 in the limit of asymptotically large energies 
if energy conservation is included in the calculation. In contrast, a two parameter 
fit of Ca/Cf and a constant offset to the gluon jet multiplicity, proposed in a recent 
study, does not yield Ca/Cf=2.25 in this limit. I apply the one parameter fit method 
to recently published data of the DELPHI experiment at LEP and determine the 
effective value of Ca/Cf from this technique, at the finite energy of the Z° boson, 
to be 1.74 ± 0.03 ± 0.10, where the first uncertainty is statistical and the second is 
systematic. 



1 Introduction 



At the basis of Quantum Chromodynamics (QCD), the gauge theory of strong in- 
teractions, are the color factors Ca and Cp, with values 3 and 4/3, respectively [I]]. 
Ca determines the relative probability for a soft gluon to couple to another gluon, 
while Cp determines the corresponding probability for a soft gluon to couple to a 
quark. The ratio Ca/Cf is perhaps the most fundamental quantity in QCD besides 
the strong interaction coupling strength, as- Currently, the most accurate measure- 
ments of Ca/Cf are from angular correlations between jets in four jet e + e~ events || 
and from the ratio of soft particle multiplicities at large transverse momenta to the 
jet axes between unbiased gluon and quark jets 0. 

Recently, a new method to measure Ca/Cf was proposed H, based on the scale 
dependence of the mean particle multiplicity in e + e~ three jet events, N 3 _j et . This 
method utilizes a next-to-leading-order (NLO[|) analytic expression for N 3 _j et ||, in 
conjunction with a constant offset term N° Q for the gluon jet multiplicity, to perform 
a two parameter fit of Ca/Cf and N°. The constant N° is intended to account for non- 
perturbative effects. The variable Ca/Cf is introduced using an analytic expression 
for the mean multiplicity ratio between gluon and quark jets, r. The expression used 
for r does not incorporate recoil effects (energy conservation). 

In this paper, I examine the determination of Ca/Cf from multiplicity in three 
jet events. My principal purpose is to test the analytic expressions for N 3 _j et and 
r. The theoretical expressions are tested by fitting them to event samples generated 
with QCD multihadronic event generators. The main conclusions are that to obtain 
the correct asymptotic result Ca/Cf=9/4=2.25 from the method it is necessary to 
use the pure QCD result [§] without the offset term N°, and to include recoil effects 
in the expression for r. As a consistency check, I apply my method to Monte Carlo 
events with Ca=Cf=4/3 to verify that the fitted result for the parameter Ca/Cf is 
consistent with unity in this case. 

Having established a fitting technique that yields the correct results in the limiting 
cases of (1) asymptotically large energies, and (2) Ca=Cf, I apply the method to 
recently published data [f| of the DELPHI experiment at the LEP Collider at CERN. 
I thereby determine the effective value of Ca/Cf from this method at the finite energy 
of the Z° boson. 



^lso referred to as MLLA. 
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2 Theoretical framework 



An analytic expression for the topology dependence of the mean particle multiplicity 
in e + e~ three jet quark-antiquark-gluon qqg events, N 3 _j et , valid in the NLO approx- 
imation of perturbation theory, is given by eq. (6.43) of ||: (see also 0): 
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where N e + e - (Q) is the mean inclusive particle multiplicity of e + e annihilation events 
at energy scale Q. The quark and gluon jet scales E* and p± are (see eqs. (6.38) and 
(6.41) of §): 
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with p q , p-g and p g the 4-momenta of the q, q and g. E* is the energy of the quark 
or antiquark in the qq rest frame, while p± is the transverse momentum of the gluon 
with respect to the qq axis in that frame. These equations are valid for massless 
quarks and gluons. 

The quantity r(Q) is the ratio of the mean multiplicities between gluon and quark 
jets. It has been calculated analytically in the next-to-next-to-next-to-leading-order 
(3NLO) approximation of perturbation theory, including recoil effects H: 



r(Q) = r (l -ri7o 



(4) 



where j (Q) = ^2CaO!s(Q) / n , r =CA/CF, and the correction terms n, r 2 and r 3 are 
constants in QCD which functionally depend on the color factors through terms pro- 
portional to 1/ro and 1/Ca M- r depends on the scale Q only through 0:5: 



(5) 



with y = ln(Q/A), A a cutoff which defines the limit of perturbative evolution, (3q = 
(IICa — 2n f )/3, (3\ = [17(Ca) 2 — n f CA(5 + 3/r )]/3, and n f the number of active quark 
flavors. In this paper, I set nf=5 and use the corresponding result for A found in a 
fit of the 3NLO expression for quark jet multiplicity || to inclusive e + e~ data. This 
result, A=0.148 GeV ||, is similar to the value of A^ |T0||.p| 



A and A-rjo are strongly related to each other but are not necessarily the same. 
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3 Analysis technique 



Three jet events are selected using standard jet finding algorithms (see section f|) • Two 
of the jets are identified as the quark (q or q) jets, the other as the gluon jet. The 
4-momentum of each jet is assigned to the underlying q, q or g. Since expressions (|2]) 
and (^|) are based on massless kinematics, the jet momenta are modified to obtain 
massless jets. First, the jets are assigned calculated energies, E ca ic, based on the 
angles between jets, assuming the jets are massless (see e.g. Second, the jet 

3-momenta are scaled as follows: 

r -> r jet— finder ; \U ) 

| P jet— finder \ 

with Pj e t-finder the jet 3-momentum determined by the jet finder. The quark and 
gluon 4-momenta defined by p=(E ca ic.,P) are used to determine the scales (|2]) and (|3]). 
This method of defining massless jets is referred to in the literature as the EO 



scheme IT2 



The values of N e + e - (2E*) and N e + e - (p±) in eq. ([[]) are determined using parametriza- 
tions of N e + e -(Q) versus Q. These parametrizations are based on sixth order polyno- 
mials in ln(Q). A parametrization is determined independently for each Monte Carlo 
event sample^ and for the data. For the Monte Carlo samples, the parametrizations 
are obtained by a fit to the predicted values of N e + e - versus Q=E CAn . in the interval 
between 10 GeV and 10 TeV, where -E c .m. is the center-of-mass (cm.) energy. For 
the data, a fit is made to measurements of N e + e - for 12 GeV < E c rxi , < 189 GeV.Q The 
parametrizations provide good representations of the multiplicity in all cases. 

The analytic expression for r (eq. (|])) is introduced into eq. ([I]). r 1; r 2 and r 3 in 
expression (f|) depend on l/r and 1/Ca, as stated above. Similarly, (3\ in ([D depends 
on l/r , while (3q and f3\ depend on Ca- The leading tq term in (||) and the l/r terms 
in ri, r 2 , r 3 and Pi form the fitted parameter. Ca in ri, r 2 , r 3 , j3 and (3\ is set equal 
to its QCD value of 3. is then determined in a one parameter fit of eq. (0) to the 
Monte Carlo or experimental results for N 3 _j et as a function of scale. The DELPHI 
Collaboration recently presented a similar study 0. I discuss the DELPHI method 
and results in section 0. 



3 Hcrwig at the partem and hadron levels, and Jetset at the parton level with Ca=Cf, see section [ 
4 The data used are the same as presented in Fig. 2 of ||. 
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4 Monte Carlo samples and event selection 



For the principal Monte Carlo based results I present, I use event samples generated 
with the Herwig Monte Carlo multihadronic event generator version 5.9. The 



parameter values used for Herwig are the same as those given in Herwig contains 
the most complete computer simulation of QCD presently available, including terms 
up to and beyond the next-to-next-to-leading order (NNLO) approximation. In this 
sense Herwig resembles an analytic calculation. In addition, Herwig implements exact 
energy-momentum conservation at each parton branching and a model for hadroniza- 
tion. Herwig yields the correct asymptotic result of 2.25 for the multiplicity ratio 
r |TJJ and related quantities H]. It provides a good description of gluon and quark 



jet properties up to the highest available e + e~ energies. Thus, Herwig provides a 
suitable QCD reference sample. Herwig generally predicts that QCD variables reach 
their asymptotic values at center of mass energies of several TeV or more, depending 
on the variable. In the following, I choose 10 TeV as the canonical cm. energy at 
which to test my fit method in the asymptotic limit. 

For studies with Ca=Cf=4/3, I employ a special version of the Jetset Monte Carlo 



muiltihadronic event generator fl5fl , version 7.4, with parameter values given in [16 
In addition to setting Ca=Cf, I turn off gluon splittings, g — > qq. The reason Jetset is 
used for these studies, and not Herwig, is that Herwig does not allow Ca=Cf- Jetset 
is based on leading order (LO) QCD with a simulation of coherence effects due to 
higher orders. The standard version of Jetset does not yield the correct asymptotic 
result for r, as seen from Fig. 2 of except perhaps at exceptionally high energies 
(£ l cm .>>100 TeV ??). Thus the QCD predictions of Jetset should be treated with 
precaution. For the present purposes it is sufficient that quark and gluon jets have 
the same internal properties, such as multiplicity, if Ca=Cf- This is satisfied by the 
special version of Jetset at the parton level. By parton level, I mean the ensemble of 
quarks and gluons which are present at the termination of the perturbative stage of 
evolution. 

Three jet events are constructed from these samples by adjusting the resolution 
scale(s) of a jet finding algorithm for each event so that exactly three jets are found. 
I choose three jet finding algorithms: the kj_ [ 17] , Jade Jl8|| and cone [l{| jet finders. 



These three algorithms are very different in their treatment of soft particles. The 
difference in the results found using the three algorithms therefore provides a rigorous 
test of the jet finder independence of the method. I retain events if the angle between 
the highest energy jet and each of the two other jets is the same to within 5°, the 
so-called "Y events. "(] An example of a Y event is shown in Fig. [TJ. Measurements of 
the particle multiplicity of Y events as a function of topology (i.e. scale) have recently 
become available 0. I wish to utilize these data for my fits (section ^). This provides 



5 Y events were first studied in 1 1 1 
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my principal motivation for selecting Y events. For Y events, the three jet event 
multiplicity N 3 _j et and the scales @ and (§) depend only on i? c . m . and one angle 
in the event, conveniently chosen to be 9\ (see Fig. [I]). For fixed E c _ m _, Q\ therefore 
determines the scale. 

For the Monte Carlo events used here, the quark and gluon jets are identified 
using parton level Monte Carlo information. The directions of the primary quark and 
antiquarkf] are determined after their perturbative evolution has terminated. The jet 
closest to the direction of the evolved primary quark or antiquark is considered to be 
a quark jet. The distinct jet closest to the other evolved primary quark or antiquark 
is considered to be the other quark jet. The remaining jet is identified as the gluon 
jet. This algorithm is applied to jets at both the parton and hadron levels. By hadron 
level, I mean the level after hadronization, with charged and neutral particles with 
lifetimes greater than 3 x 10~ 10 s treated as stable. Hence charged particles from 
the decays of K° and weakly decaying hyperons are included in the definition of the 
hadron level multiplicity. 



5 Monte Carlo based results 



I begin by studying Herwig events at the parton level, with E c _ m _=10 TeV. This 
large energy is meant to ensure that the fit results are asymptotic, as mentioned 
above. The mean multiplicity of these events as a function of the opening angle Q\ 
is shown in Fig. ^|a. The results are shown for the three jet algorithms. The results 
of the three algorithms are seen to be similar for angles larger than about 80°. As 9\ 
becomes smaller, the two lower energy jets are not as well separated and background 
from two jet-like events increases. Different jet finders have different efficiencies for 
selecting background: thus the results of the jet finders diverge. Since the results 
should be independent of the choice of a jet algorithm to be sensible, I restrict the 
fits to the region where the predictions of the jet finders approximately agree, namely 
80°<fli<120°. 

The solid curve in Fig. 0a shows the result of the one parameter fit of eq. (|l|) to the 
event multiplicity determined using the k± jet finder. The curve provides a reasonable 
description of the multiplicity inside the fit region. Outside this region, i.e. for angles 
less than 80°, the fitted curve does not describe the event multiplicity well. The 
multiplicity of the events in Fig. 0a is not well defined for 9i < 80°, however, since 
the results from different jet finders disagree strongly. Therefore I do not consider the 
discrepancies between the fitted curve and the jet finder based results for 9 1 < 80° to 

6 i.e., the q and q produced directly in the electroweak decay of the virtual Z°/7 in e + e~ — > 
(Z°/7)* — ► hadrons events. 
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be meaningful. 



The results for r are summarized in the top portion of Table |l[ Taking the result 
found using the k± jet finder as the central value, with half the difference between 
the extreme values found using the different jet finders as a systematic uncertainty, 
I obtain r =2.248 ± 0.010 (stat.)±0.024 (syst.), consistent with the QCD asymptotic 
expectation of 2.25. 

The analogous results for Jetset at the parton level with Ca=Cf are shown in 
Fig. ^d. Again, the cm. energy is 10 TeV. The predictions of the three jet finders are 
seen to be similar only for $i > 90°. Therefore, I limit the fit range to 90° < 0\ < 120° 
in this case. The results for r$ are given in the bottom portion of Table |I|. Combining 
the results in the manner described in the previous paragraph yields ro=1.012 ± 
0.009 (stat.)±0. 027 (syst.), consistent with unity. 

Thus a one parameter fit of tq yields the correct results in the limiting cases of 
(1) QCD at asymptotically large energies and (2) Ca=Cf, as long as the fit range is 
restricted to regions where the results of the different jet finders agree, or, equivalently, 
to regions where the fitted curves provide a good description of the multiplicity. The 
fits generally yield y 2 /d.o.f. ~ 1 (Table p. 

It is of interest to determine the importance of energy conservation in the ex- 
pression for r. To this effect, I replace eq. (|J) by the corresponding result in the 
NNLO approximation both with |2(| and without || recoil effects, and repeat the fit 



of Herwig events described above. The NNLO approximation is used for this test, 
and not the 3NLO approximation, because a 3NLO expression for r without energy 
conservation is not available. For the k± based event sample, the NNLO calculation 
with energy conservation yields ro=2.254 ± 0.010 (stat.), not very different from the 
3NLO result presented above. The corresponding result without energy conservation 
is only 2.079 ± 0.009 (stat.), however, significantly smaller than 2.25. This implies 
that it is important to include energy conservation in the analytic expressions, even 
for E c , m =10 TeV. 



In Fig. |3|, I show the fitted results for r as a function of E cra , using Herwig events 
at the parton and hadron levels. The events are selected using the k± jet finder. 
The hadron level multiplicity is based on charged particles only. The fit interval is 
80° < 6i < 120°, i.e. the same as in Fig. |2|a. The fitted curves provide good descriptions 
of the multiplicity within the fit region for all cm. energies at both the parton and 
hadron levels. The fit results are observed to have only a moderate dependence on the 
choice of the jet algorithm, generally similar to that indicated in Table [I] for parton 
level events or in item 1 of the list presented below in section |6] for hadron level events. 
From the parton level curve (solid line) it is seen that the asymptotic result tq ~ 2.25 
is reached for -E c . m . ~3 TeV. The hadron level curve (dashed line) converges to the 
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asymptotic limit much more slowly, however. As a consequence, the hadronization 
correction, defined by the ratio of the parton to the hadron level results, is fairly 
large. The hadronization correction is predicted to be 1.30 at the mass of the Z° and 
1.06 at 10 TeV. 

The dotted curve in Fig. ^ shows the fitted results for tq at the parton level if the 
NNLO expression for r without recoil effects is used in place of the 3NLO expression. 
The QCD asymptotic limit of 2.25 it not attained in this case, again emphasizing the 
importance of energy conservation. 



6 Data based results 

Recently, the DELPHI experiment at LEP presented measurements of the charged 
particle multiplicity of Y events and the scales (fj) and (|3|) as a function of the opening 
angle 6i [|J. The results are based on the k± jet finder with E cm =91 GeV. These 
data allow me the possibility to determine the effective value of r at the scale of 
the Z° using my one parameter fit method. The DELPHI multiplicity measurements 
are shown in Fig. |4j. The result of the one parameter fit is shown by the solid curve. 
The fit range employed is 78° <6 1 < 120°, similar to the interval of 80° < Q x < 120° 
I choose for Herwig events (Figs. [|a and |3|). The small difference in the choice of 
fit interval between the Herwig and DELPHI samples is not important (see item 2 
below). The analytic curve provides a good description of the measurements within 
the fit region, yielding \ 2 /d.o.f. =8.9/8. 

The result for the fitted parameter is ro=1.737 ± 0.032 (stat.). To estimate a 
systematic uncertainty for this result, I consider the following. 

1. Jet finder dependence. The DELPHI results are presented for the k± jet finder, 
but not for the Jade or cone jet finders. Herwig at the hadron level with 
E c . m =91 GeV yields r =1.585 for the k± jet finder, 1.601 for the Jade jet 
finder, and 1.516 for the cone jet finder, where the statistical uncertainty is 
0.008 in all cases. Half the difference between the extreme values is taken as a 
systematic uncertainty. 

2. Fit interval. The fit interval I choose for the DELPHI data is 78° < 6>i < 120°, as 
stated above. Decreasing the lower limit of this interval to 60° yields ro=1.705± 
0.025 (stat.), while decreasing the upper limit to 90°, with the lower limit at the 
standard value, yields r =1.755±0.038 (stat.). I take half the difference between 
these values as a systematic uncertainty. A further test of the choice of the fit 
interval is presented in Fig. |5]. Fig. |5|a shows the fitted results for 
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function of 6f n -, where 6f n - is the lower limit of the fit range 6f n - <6i< 120°. 
The corresponding values of \ 2 /d.o.f. are shown in Fig. [|b. The \ 2 /d-o.f. is 1.4 
or less for Of 1111 ' > 60° but much larger for 9™ m ' < 60°. This provides justification 
for not extending the fit range below 60°, i.e. the fit is restricted to an interval 
where the theoretical expression eq. (|l|) describes the data accurately. 

Parametrization o/N e + e - versus Q. Rather than use a polynomial parametriza- 
tion of N e + e - versus Q (section D), I use the parametrization based on the 3NLO 
expression for quark jet multiplicity || with the parameter values in |9[]. This 
yields r =1.804 ± 0.032 (stat.). The difference with respect to the standard re- 
sult is taken as a systematic uncertainty. Note that the polynomial provides a 
better description of N e + e - versus Q than the 3NLO expression. 



4. Value of A. Setting A in eq. (|) to the PDG value of A M g=0.220 GeV flO 
rather than using 0.148 GeV (section ^), yields r =1.761 ± 0.033 (stat.). I take 
the difference with respect to the standard result as a systematic uncertainty. 

5. Averaging procedure for E* and p±. The DELPHI results for the quark and 
gluon scales eqs. @ and (|[) are found by taking the geometric means of E* 
and p±, averaged over the event sample, as a function of 9\. For the Monte 
Carlo based results presented in section |5], I employ the much more common 
arithmetic means. For hadron level events at 91 GeV, Herwig with the k± jet 
finder yields r =1.629 for geometric means and r =1.585 for arithmetic means, 
where the statistical uncertainty is 0.008 in both cases. The difference between 
these values is taken systematic uncertainty. 

6. Number of active flavors, nf. Using nf=3 rather than nf=5 in the analytic 
expressions for r and as (eqs. (f|) and (Q)), and correspondingly evaluating as 
using A=0.322 GeV @ rather than A=0.148 GeV, yields r =1.735±0.029 (stat.). 
The difference with respect to the standard result is taken as a systematic 
uncertainty. 



The systematic uncertainties are summarized in Table |2|. The largest terms arise from 
the parametrization of N e + e -, the averaging procedure for the scales, and the choice 
of jet finder, in that order. The terms are added in quadrature to define the total 
systematic uncertainty. The final result for the effective value of ?"o=Ca/Cf at the 
scale of the Z° is: 

r = 1.737 ± 0.032 (stat.) ± 0.097 (syst.) . (7) 

Multiplying this value by the hadronization correction of 1.30 mentioned at the end 
of section | yields r = 2.26 ± 0.04 (stat.) ± 0.12 (syst.) as the corresponding result at 
the parton level, consistent with the QCD asymptotic prediction of 2.25. The data 
based results I obtain at the hadron and parton levels are shown by the solid and 
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open points in Fig. 0. The experimental results lie somewhat above the Herwig curves 
but are generally consistent with them. 



7 Comparison to a two parameter fit method 



In their recent publication 0], the DELPHI Collaboration presented an alternative 
method to determine to—Ca/Cf from three jet event particle multiplicity. I used the 
data of that study, shown in Fig. |], to obtain the results of section |[ The DELPHI 
analysis is based on a fit of the expression 



N 3 _ jet = N e+e -(2E*)+r(p J 



N e + e -(p±) n q" 

9 S 



to the three jet event multiplicity data, where N° is a parameter meant to account for 
differences in the hadronization of gluons and quarks. The DELPHI analysis differs 
from mine principally by using expression @ rather than expression ([!]), by using the 
NNLO result for r without energy conservation || rather than the 3NLO expression, 
and by invoking a two parameter fit of r and N° rather than a one parameter fit. 
The DELPHI analysis also differs from mine in the choice of fit range and in the 
parametrization of N e + e - versus Q. In the DELPHI study, the fit range is 30° < Q\ < 
120° rather than 78° < Q\ < 120° and the parametrization of N e + e - versus Q is based 
on the NLO expression for quark jet multiplicity in e + e~ annihilations Klj rather 
than a polynomial. The DELPHI results utilize events at E cm =91 GeV selected 
using the k± jet finder, as stated in section ||. 

Repeating the DELPHI analysis, viz. a two parameter fit of eq. (||) to the data in 
Fig. |], using the DELPHI values of E* and p±, the expression for r in 0, a fit range 
from 30° to 120°, and the NLO expression for quark jet multiplicity to parametrize 
N e + e - versus Qf\ I obtain 

r = 2.200 ± 0.066 (stat.) ; (9) 
Nj = 1.46 ±0.10 (stat.) . (10) 

The x 2 /d.o.f. of the fit is 13.2/16. The results (^) and (|10D are very similar to those 
of DELPHI, namely r =2.251± 0.063 (stat.) and 1.40 ±0.10 (stat.) @. The result (|) 
for r is shown by the solid point in Fig. ^. 

The value of r derived from the DELPHI two parameter fit method is numerically 
very similar to the QCD asymptotic result Ca/Cf=2.25. On this basis, DELPHI 
suggests M that their analysis provides a measurement of that quantity. To test 



7 For the NLO parametrization of quark jet multiplicity, I use the parameters in p2| 
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this hypothesis, I determine the results of the DELPHI method in the two limiting 
cases discussed in section [|: (1) asymptotically large energies and (2) Ca=Cf- For 
Herwig events at the parton level with E crai =10 TeV, the DELPHI fit method yields 
r =2.80±0.03 , 3.19±0.03 and 4.03 ±0.05 for events selected using the k±, Jade and 
cone jet finders, respectively, where the uncertainties are statistical. These values are 
much larger than 2.25 and exhibit a strong dependence on the choice of the jet finder, 
in contrast to the results of section |5] (top portion of Table [I]). The analogous results 
for Jetset at the parton level with E cia =10 TeV and Ca=Cf are 1.76 ± 0.03 (stat.), 
2.05 ±0.02 (stat.), and 2.68 ±0.03 (stat.), which are inconsistent with unity and again 
exhibit a strong jet finder dependence. This is also in contrast to the results of 
section [| (bottom portion of Table [I]). For the above results, the NLO expression for 
quark jet multiplicity is fitted to the MC predictions of N e + e - versus £ c . m . for both 
the Herwig and Jetset samples, using scale values between 20 GeV and 10 TeV. The 
results are similar if the polynomial parametrizations discussed in section ^ are used 
instead. 

The dashed and solid curves in Fig. ^| show the results I obtain for r from applying 
the DELPHI fit method to Herwig events at the hadron and parton levels. The 
results are shown as a function of E cm . The hadron level results are based on 
charged particles only. The event samples are selected using the k± jet finder. Thus 
Fig. P is the analogue for the DELPHI method of the results I show in Fig. || for my 
method. The hadron level curve in Fig. |6| is seen to be in general agreement with 
the experimental result @ at the scale of the Z°. Asymptotically, the hadron level 
prediction reaches a value of about 2.7, however, much larger than 2.25. The parton 
level curve exceeds 2.25 by a large margin even at E c , m =91 GeV. 

On the basis of the above results, I conclude that the DELPHI fit method probably 
does not measure Ca/Cf, and that the similarity of the hadron level result (|9]) to 
the asymptotic prediction Ca/Cf=2.25 is most likely a coincidence. As a last note I 
remark that if energy conservation is included in the NNLO expression for r, the result 
(P) increases to ro=2.479 ± 0.081 (stat.). Thus if energy conservation is incorporated 
into the DELPHI fit method, the value of vq derived from charged hadrons at 91 GeV 
is no longer similar to 2.25. 



8 Summary and Conclusions 

In this paper, I have presented a test of the QCD expression for the topology (scale) 
dependence of particle multiplicity in e + e~ three jet events. Using event samples 
generated with the Herwig Monte Carlo multihadronic event generator as a reference, 
I find that the QCD expression yields the correct result Ca/ Cp=2.25 in the asymptotic 
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limit of large energy scales Q ~ 3 TeV as long as it is used in conjunction with an 
expression for r which incorporates energy conservation, where r is the ratio of mean 
particle multiplicities between gluon and quark jets. This emphasizes the importance 
of energy conservation in QCD analytic expressions, even at large scales. My analysis 
is based on a one parameter fit of Ca/Cf to three jet event mean particle multiplicity 
function of the topology of the event. 

As a second test, I apply my method to a sample of Monte Carlo three jet events 
in which the color factors are set equal, Ca=Cf- I obtain Ca/Cf ~1 in this case, 
demonstrating the self consistency of the technique. 

Applying my fit method to recently published data [[| of the DELPHI experiment 
at LEP, I obtain C A /C F =1.737 ± 0.032 (stat.) ± 0.097 (syst.) as the effective value 
of the color factor ratio at E c . m =9l GeV from this technique. This result is based 
on charged particles. It is of interest to compare this result to related measurements 
at E cra =91 GeV based on the charged particle multiplicity ratio between gluon 
and quark jets, r c ^.. The experimental result for r c ^. in full phase space is 1.514 ± 
0.019 (stat.) ± 0.034 (syst.) |H| . The corresponding result for r ch . in limited phase 
space, defined by soft particles with large transverse momenta to the jet axes, is 
2.29 ± 0.09 (stat.) ± 0.15 (syst.) P|. All these measurements - the one presented here 
and the two based on r c h. - are predicted to equal 2.25 in the limit of large energies. 
The result presented here is seen to be intermediate to the two based on r ch , both 
in value and in the size of the uncertainty. Whereas r c h. in limited phase space has 
already attained its asymptotic value at E c , m =91 GeV, r c h. in full phase space and 
the result presented here are sub- asymptotic at this scale. 

Last, I test a two parameter fit method to determine Ca/Cf from particle mul- 
tiplicity in three jet events, proposed in a recent study ||. I find that this method 
does not yield the correct results Ca/Cf ~ 2.25 or Ca/Cf ~1 in the two limiting 
cases of QCD at asymptotic energies or identical color factors Ca=Cf, in contrast to 
my method. Thus, I conclude that this two parameter fit method probably does not 
measure the color factor ratio. 
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r =C A /C F xVd.o.f. 


Herwig partons 
kj_ jet finder 

TpHp lpf fiTiHpr 

Cone jet finder 


2.248 ±0.010 5.9/8 
2 269 ± 010 12 8/8 
2.221 ±0.013 14.7/8 


Jetset partons, Ca=Cf 
kj_ jet finder 
Jade jet finder 
Cone jet finder 


1.012 ±0.009 5.9/6 
1.032 ±0.007 3.9/6 
0.979 ±0.008 6.7/6 



Table 1: Results of a one parameter fit of r =CA/CF to the parton level multiplicity 
in three jet events, as predicted by the Herwig QCD multihadronic event generator, 
and by the Jetset multihadronic event generator with Ca=Cf- The -E c .m. value is 
10 TeV for both samples. The fits are performed using the 3NLO expression for the 
multiplicity ratio between gluon and quark jets, r. The uncertainties are statistical. 



Systematic term 


Ar 


1. Jet finder dependence 


0.043 


2. Fit interval 


0.025 


3. Parametrization of N e + e - 


0.067 


4. Value of A 


0.024 


5. Averaging procedure for E* and p± 


0.044 


6. Number of active flavors 


0.002 


Total 


0.097 



Table 2: Summary of systematic uncertainties for the effective value of tq at the scale 
of the Z° as determined using data [|J. 
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Figure 1: Schematic representation of a three jet qqg event produced in e + e~ anni- 
hilations with a Y event topology |TTJ , in which the angle between the highest energy 
jet and each of the two lower energy jets is about the same. The angle Q\ opposite 
the highest energy jet is used to specify the event topology. 
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( b ) Jetset partons, C A =C F ,E c m =10 TeV 
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Figure 2: (a) The mean parton level multiplicity of three jet Y events as a function of 
the opening angle 9\, for events generated with the Herwig multihadronic event gen- 
erator, (b) The analogous results for events generated with the Jetset multihadronic 
event generator with Ca=C F - The event samples in (a) and (b) are selected using 
the k±, Jade and cone jet finders. The cm. energy is 10 TeV. The solid curves show 
the results of a one parameter QCD fit to events selected using the k± jet finder. The 
fits are performed within the regions shown. 



16 



2.8 



2.6 



2.4 - 



2.2 



U 2 h 



1.8 - 



1.6 - 



1.4 



1.2 - 




/ 



/ 



Data, hadron level 

Data, parton level 

Herwig, hadron level 

Herwig, parton level 

Herwig, parton level, 
no recoil effects in r 



10 



E. 



cm. 



10 

(GeV) 



10 



Figure 3: Results of a one parameter fit of r =CA/CF as a function of the cm. 
energy, for Herwig Monte Carlo three jet Y events at the parton and hadron levels. 
The events are selected using the k±_ jet finder. The corresponding results for data [|] 
at E cm =91 GeV are shown by the open and solid points. The hadron level results 
are based on charged particles only. For the data points, the vertical lines show the 
total uncertainties, with statistical and systematic terms added in quadrature. The 
small horizontal lines indicate the statistical uncertainties. 
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Figure 4: Measurements ||J] of the mean charged particle multiplicity of three jet 
Y events as a function of the opening angle 6\, for E cm .=91 GeV. The events are 
selected using the k± jet finder. The solid curve shows the result of a one parameter 
fit to the data within the fit region shown. 
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Figure 5: (a) Results of a one parameter fit of r =CA/CF to measurements (4| of the 
charged particle multiplicity of three jet Y events at E c , m =91 GeV, as a function of 
the lower limit #™ m ' of the fit range Of 1111 ' <6\< 120°. (b) The corresonding values of 
X 2 /d.o.f. The events are selected using the k±_ jet finder. 
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Figure 6: Results for ro=CA/Cp from a two parameter fit method Q as a function of 
the cm. energy, for Herwig Monte Carlo three jet Y events at the parton and hadron 
levels. The events are selected using the k±_ jet finder. The corresponding hadron 
level result found using data j| at i£ cm .=91 GeV is shown by the solid point. The 
hadron level results are based on charged particles only. The uncertainty shown for 
the data point is statistical. 
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